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Perovskite Manganites Hosting Versatile Multiferroic Phases 
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Complete magnetoelectric (ME) phase diagrams of orthorhombic _RMn03 with and without mag- 
netic moments on the R ions have been established. Three kinds of multiferroic ground states, the 
o6-cycloidal, the fec-cycloidal, and the collinear _E-type phases, have been identified by the distinct 
ME responses. The electric polarization of the _E-type phase dominated by the symmetric spin 
exchange (Si ■ Sj) is more than 10 times as large as that of the bc-cycloidal phase dominated by the 
antisymmetric one (Si x Sj), and the ME response is enhanced near the bicritical phase boundary 
between these multiferroic phases of different origins. These findings will provide an important clue 
for the development of the magnetically induced multiferroics. 

PACS numbers: 75.30.Kz, 77.80.-e, 75.80.+q 



It has been a long standing problem how to enhance 
the correlation between magnetism and ferroelectricity in 
a solid. Possible solutions can be found in recent stud- 
ies on the magnetically induced multiferroics PHH, which 
can be classified into two types; one is driven by antisym- 
metric exchange striction in the cycloidal spin structure, 
typified by TbMn0 3 @, H Q, and the other is driven 
by symmetric exchange striction in the commensurate 
collinear spin structure. Since orthorhombic (o-) per- 
ovskite i?Mn03 (R = rare earth and Y) contains the 
both types, it provides an ideal laboratory to compare 
the respective magnetoelectric (ME) properties and ex- 
tract the essential ingredients for the development of the 
magnetically induced multiferroics. 

The ferroelectricity with a cycloidal spin order in o- 
i?Mn0 3 has been found for R = Gd, Tb, Dy, and 
Eui-^Ya; 0, 0, 0, Q. Upon the application of mag- 
netic field B along the a axis in the afr-cycloidal phase 
(in Pbnm notation) of Eui_ I Y a; Mn03, the spin-cycloidal 
plane rotates from ab to be, accompanied by the po- 
larization P rotation from a to c [3, The ori- 
gin of the ferroelectricity has been discussed in terms 
of the spin-current model or the inverse Dzyaloshinskii- 
Moriya (DM) interaction represented by the relation, 
P ~ Y, Ae i] x {Si x Sj) [TO. 03- 03. in which e tJ is 
the unit vector connecting the neighboring spins (Si and 
Sj) and both the spin-orbit and super-exchange interac- 
tion are relevant to the coefficient A. On the other hand, 
o-i?MnC>3 with R = Ho, Tm, Yb, Lu shows a commen- 
surate collinear spin order with a propagation vector q — 
(0, 1/2, 1), which is so-called E-type antiferromagnetic 
order 0, M, ES 03 • This phase is allowed to possess 
P along a due to symmetric exchange striction, which 
is independent of the spin-orbit interaction. Sergienko et 
al. proposed the emergence of substantially large P up to 



0.12 C/m 2 in the E-type phase with considering the fer- 
romagnetic nearest-neighbor interaction J\ mediated by 
e electrons as a major source of the exchange striction 



18, 1 



The rich variety of magnetic phases of 0-RM11O3 re- 
flects the J\-J% competition that varies as a function 
of i?-ion radius under the staggered orbital ordering of 



3x 2 



and 3y 



type orbitals; As GdFeGvtype lat- 



tice distortion becomes large, the antiferromagnetic next- 
nearest-neighbor interaction J2 in the ab plane becomes 
competitive with the ferromagnetic interaction J\. So 
far, ME phase diagrams of 0-RM11O3 have been made for 
R from La to Dy or EU0.5Y0.5 (see Figs. 1(a) and 1(b)) 
[SHIMI!]],^. However, since o-i?Mn0 3 with a smaller 
R ion than Dy needs a high pressure (HP) technique to 
synthesize, a complete ME phase diagram of o--RMnC>3 
including the neighboring area between the cycloidal and 
the E-type phases is absent. In fact, although the fer- 
roelectricity in the E-type phase has been confirmed for 
o-i?Mn0 3 with R = Y, Ho m 



and Tm [17(, intrinsic ME 
properties inherent to the Mn-spin arrangement alone re- 
main unclear because of the intervention by the magnetic 
R ions. Besides, the reported P values vary widely de- 
pending on the materials. As a result, quantitative es- 
timations of P of the E-type phase and its microscopic 
origin are under intensive debate. 



This Letter reports complete ME phase diagrams of 
a series of o-i?MnC>3 with nonmagnetic R ions (R = 
Eui-^Ya; and Yi^Lu^) together with the system con- 
taining the magnetic R ion (R = Dy, Ho, Er, Tm, 
Yb). By using high-quality polycrystalline samples of 
o-i?Mn03, we have confirmed substantially large P of 
nearly 5000 /iC/m 2 for the E-type phases, yet one order 
of magnitude smaller than the predicted values lH, li| , 
and clearly demonstrated the generic transition of the 
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FIG. 1: (Color online) Phase diagrams of o-RMnOs with (a) 
magnetic R = Gd, Tb, Dy, Ho, Er, Tm, Yb, and (b) non- 
magnetic R = Eui-jYx and Yi-^Luj, (PM, IC, and AF de- 
note paramagnetic, incommensurate, and antiferromagnetic, 
respectively), (c) polarization P value at 2 K for the com- 
pounds with nonmagnetic 7? ions (P HP , which is multiplied 
by a calibration factor (= 6) [26|], is for P of the polycrystals, 
and P a is for P along a of the single crystals taken from ref. 
[jj, (d) k in the magnetic pro pagation vector q = (0, k, 1) 
(taken from refs. [li, Q3, [H, EE 03, [H, El), and (e) unit 
cell volume as a function of the 7?-ion radius, of which coor- 
dination number is assumed to be 8 for an orthorhombically 
distorted perovskite lattice. Tni was determined by magnetic 
and dielectric measurements. Tn2 and Tn3 were determined 
from the measurements of P. Tni and Tn2 for i?Mn03 with 
R = Eu, Gd, and Tb were taken from ref. 0]. The shaded 
area in (b) represents a possible phase-coexisting region. 



multiferroic ground state from the afo-cycloidal, to bc- 
cycloidal, and eventually to the E-type phase upon de- 
creasing the size of the R ion. Furthermore, we found an 
enhanced ME response characteristic of a bicritical phase 
boundary formed by the different multiferroic phases. 

Polycrystalline samples of o-i?Mn03 with R = Dy and 
R not larger than Y were synthesized under HP by us- 



ing hexagonal (h-) i?MnC>3 as precursors. First, single 
crystals of h-RMnOz were prepared by a floating-zone 
method, and then heat-treated for 1 hr in the range of 
1323 - 1373 K under a HP of 5.5 GPa. Just for com- 
parison between the single crystalline and polycrystalline 
specimens, polycrystalline samples of o-Eui-^Y^MnOs 
with x = 0.2, 0.4, 0.6 were synthesized by heat treat- 
ment under HP using polycrystalline precursors prepared 
by grinding the single crystals of the o-phase. The grain- 
boundary effects on ME measurements are expected to 
be minimized owing to the high-quality precursors. How- 
ever, as for o-Eui-zY-jMnOa with x = 0.75, 0.9, the at- 
tempts to grow crystals gave the mixtures of the h- and 
the o-phases, which were adopted as the precursors. The 
linear change in the unit cell volume as a function of 
the i?-ion radius ensures the successful syntheses of a se- 
ries of o-RMnOs (Fig. 1(e)). For measurements of P, 
gold electrodes were deposited on the polished faces (4 
mm 2 ) of the platelet samples with typically 0.2 mm in 
thickness. As a poling procedure, an electric field of 800 
V/mm was applied at 40 K, followed by cooling to 2 K. 
The displacement current was measured with increasing 
temperature at a rate of 5 K/min or sweeping B at a rate 
of 100 Oe/sec, and was integrated as a function of time 
to obtain P. Magnetic susceptibility was measured by a 
SQUID magnetometer. 

Figure 2 shows temperature dependence of P of o- 
i?MnC>3 prepared by the HP technique, with high- 
lighting different responses of P to external B for o- 
Euo.iY .gMn03 and o-LuMn03. As exemplified in Figs. 
2(a) and 2(b), the compounds undergo a transition to 
an incommensurate (IC) sinusoidal phase at Tni where 
X takes a maximum. Then, they show a second tran- 
sition to the .E-type phase at Tn2 where P sets in. In 
Figs. 2(a) and 2(b), P in zero magnetic field is plot- 
ted against temperature for nonmagnetic and magnetic R 
ions, respectively. The compounds with R = Eun.iYn.g, 
Yi-yLuy, Ho, Er, Tm, Yb are supposed to possess the 
E-type phase, and all the compounds but for R = Er and 
Eun.iYn.g exhibit fairly large P of about 800 nC/m 2 at 
the lowest temperature. By contrast, R = Evli- x Y x [x 
= 0.4, 0.6, 0.75) and Dy sample possess relatively small 
P values, implying the different mechanism of P genera- 
tion. We should note here that o-YMn03 and o-ErMn03 
show the stepwise temperature dependence of P. The 
neutron diffraction studies have suggested the IC q vec- 
tors of (0, 0.435, 1) and (0, 0.433, 1) for the magnetic 
ground states in o-YMn03 and o-ErMn03, respectively 



2J, |25j. Thus, these stepwise changes may indicate the 
presence of a small amount of the cycloidal or the sinu- 
soidal phase which competes and thus coexists with the 
E-type phase. The signature for such a phase coexistence 
is explicitly found for o-Eun.iYrj.gMnOs showing smaller 
P and the larger effect of applied B on P in compari- 
son with LuMnOs (see Figs. 2(c) and 2(d)). The possi- 
ble phase coexistence in o-YMn03 and 0-Eun.iY n .gMn03 
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FIG. 2: (Color online) Temperature dependence of polar- 
ization P of o-J?Mn03 with (a) nonmagnetic R (Eui-jYj, 
Yi-yLuj,) and (b) magnetic R (Dy, Ho, Er, Tm, and Yb) 
in the absence of external magnetic field B, and (c) P of o- 
Euo.iYo. 9 Mn03 (multiplied by 2) and (d) P of o-LuMn0 3 
measured in external B perpendicular to the applied electric 
field and x measured in external B of 0.1 T on increasing 
temperature. B was applied during the poling procedure as 
well as the measuring process. 



can be attributed to proximity to the first order phase 
boundary, as discussed later. As for o-ErMnC>3, a long 
range order of the E-type phase is likely to be amenable 
to the large magnetic moments of Er ions with planar 
anisotropy, leading to the stepwise change in P. 

Based on the polarization and magnetization measure- 
ments, we established the ME phase diagrams for a series 
of o-i?Mn03 with magnetic and nonmagnetic R ions, re- 
spectively (Figs. 1(a) and 1(b)). Here, one may find 
that the variation of Tni as a function of the ionic ra- 
dius of R is well systematic when the R ions are non- 
magnetic, whereas it is rather irregular when the R ions 
are magnetic. This difference suggests that the presence 
of the magnetic moment of the R ions gives substan- 
tial influence on the transition temperatures as well as 
the relative stability of the two kinds of cycloidal phases 
and the -E-type phase. Figure 1(c) shows a plot of P 
at 2 K multiplied by a correction factor (= 6) for the 
compounds with nonmagnetic R ions t 26], together with 
P a in Eui_ I Y a; Mn03, as a function of the i?-ion radius. 
Contour plot of the calibrated value of P is shown in 
Fig. 3(c). In the E-type phase, P reaches nearly 5000 
/iC/m 2 , which is more than 10 times as large as P (= 
384 yitC/m 2 ) in Euo.25Yo.7sMn03, and almost indepen- 
dent of the i?-ion radius except for Euo.iYo.gMn03. In 
the cycloidal phases, on the other hand, P is critically 
dependent on the i?-ion radius. According to the inverse 
DM model, P is proportional to sin9, where 9 denotes the 
angle between the nearest neighbor Mn spins. When the 
afr-cycloidal spins are of ideal cycloid and k is less than 
0.5 (9 < 7r/2), P (oc sin0) is expected to increase mono- 
tonically, as indicated by the broken line in Fig. 1(c). 
However, in the range of 0.4 < x < 0.75, P decreases as 



increasing x. This is probably because the ellipticity of 
the cycloidal spins becomes larger 22 1, or a fraction of 
the frc-cycloidal phase with smaller P becomes larger. By 
further increasing x, P shows a local minimum at x = 
0.75, where the frc-cycloidal state is the magnetic ground 
state, followed by an abrupt increase that indicates the 
emergence of the E-type phase. 

In order to reveal the difference in ME response of 
each phase of o-RMnO^ with nonmagnetic R, we have 
investigated B dependence of P at representative points 
(T = 2 K) marked with open circles in Fig. 3(a). In 
Eu .iY . 9 MnO3 and LuMn0 3 (Figs. 3(d) and 3(e)), the 
normalized P shows nearly quadratic behavior as a func- 
tion of B and little dependence on the B direction with 
respect to the P direction. The results are in accord with 
the E-type phase as the magnetic ground state for both 
compounds, because P of this phase should be propor- 
tional to Si ■ Sj and thus to 1-aB 2 when the magnetiza- 
tion is proportional to applied B. As clearly seen in Fig. 
3(a), the R- ion dependence of the transition tempera- 
ture Tn2, below which the system becomes ferroelectric, 
makes a V-shaped feature near the boundary between the 
&c-cycloidal and the E-type phases. This feature reminds 
us of the bicriticality of two competing phases separated 
by the first-order phase boundary in the presence of weak 
randomness |27(] . The 6c-cycloidal and the E-type phases 
are likely separated by the first-order phase boundary 
[28j |. and Eui_ x Y x Mn03 has also weak randomness in- 
troduced via alloying the A site. Since Euo.iY .gMn03 
is located on the verge of the transition from the E-type 
phase to the frc-cycloidal phase, the system should be sub- 
ject to the phase coexistence, leading to the enhancement 
of the effect of B on P (Fig. 3(d)) and the suppression 
of P (Fig. 2(c)). As a hallmark of the phase coexistence, 
the large hysteresis was observed in the P-B curve of 
Euo.iY .gMn03, suggesting that the decrease of P is not 
only due to the cant of spins in the E-type phase but 
partly due to the increase in the volume fraction of the 
cycloidal phase which is dominated by the Zeeman en- 
ergy gain (oc x^ 2 )- In accord with this expectation, x 
is larger for Euo.iYo.gMnOs containing a small amount 
of the large- x cycloidal phase than for LuMnOs with the 
purely E-type order (see Figs. 2(c) and 2(d)). 

On the other hand, in the <x6-cycloidal phase of 
Euo.6Yo.4MnC-3 (Fig. 3(b)) and the frc-cycloidal phase 
of Euo.25Yo.75MnC"3 (Fig. 3(c)), P changes drastically 
at critical B of about 5 T and 3 T, respectively, which 
signifies the occurrence of the P rotation. In fact, the 
study of a single crystal of Euo.6Yo.4Mn03 showed the 
5-induced P rotation from P a to P c at 4.5 T (B || a) [§]. 
However, because our samples are polycrystals, we have 
to consider the relative orientations between B, P, and 
the crystallographic axes of a single domain. To explain 
the drastic change in P under B, we considered specific 
situations which satisfy the conditions for the occurrence 
of the 5-induced P rotation (hereafter we formally use 
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FIG. 3: (Color online) (a) Corrected magnitude of polar- 
ization P displayed as a contour-plot in the phase dia- 
gram of o-RMnOa with nonmagnetic 7? (Eui-^Yx, Yi_ H Lu H ). 
Magnetic-field B dependence of normalized P at 2 K for (b) 
Euo.6Y .4Mn0 3 , (c) Eun^sYn.ysMnOs, (d) Euo.iYn.gMnOa, 
and (e) LuMnOa (P was normalized by the value at 2 K). 
Schematic illustrations of the spin rotation from the be cy- 
cloidal to the ab cycloidal (f) under B \\ c, (g) under B \\ P m 
deviating from c to a by <f>, and (h) under B _L P m deviating 
from c to a by <j> (<j> < 7r/4). These illustrations depict the 
spin rotation in a single domain both for lattice and spin. The 
relational expressions between P m or P' m and (j> are shown. 



P m and P' m for measured P before and after the P rota- 
tion, respectively). The panels in Fig. 3(f) describe the 
P rotation from P c to P a upon the application of B along 
c. After the P rotation, two kinds of domains with pos- 
itive and negative P a coexist in an equal amount (only 
for the P rotation to positive P a is shown in Fig. 3(f)), 
so that P' m should be zero. When B is slightly deviated 
from c to a by <p, finite P' m remains after the P rota- 
tion as shown in Figs. 3(g) and 3(h), because the plane 
of spin spiral inclines to be perpendicular to B to gain 
the Zeeman energy during the P rotation (see the mid- 
dle panels), and thereby one of the rotation directions is 
selected, as reported in Euo.55Y .45Mn03 [13]. If B is 
slightly deviated from c to b, P' m should be zero similarly 
to the case in Fig. 3(f). Here, the normalized P after 
the P rotation (P^/P m ) for B \\ P m and B _L P m can 
be expressed as P' m jP m = -t&n<j)P a /P c and P' m jP m = 



tan(7r/2 — cj))P a /P c , respectively. Taking the inequations 
P a > P c and <f> < 7r/4 into account, the drastic decrease 
(and even sign reversal) of P for B \\ P m and the in- 
crease of P for B ± P m shown in Fig. 3(c) can be well 
explained. The same explanation is applicable for the P 
rotation from P a to P c (Fig. 3(b)). 

In summary, we have revealed the magnetoelectric 
properties of a whole series of multiferroic o-i?Mn03 with 
both nonmagnetic and magnetic R and estimated the 
genuine values of P (~ 5000 /iC/m 2 ) in the -E-type phase, 
which is more than 10 times as large as that of the bc- 
cycloidal phase, yet one order of magnitude smaller than 
the predicted values li], 1|| ■ Furthermore, we found the 
bicritical feature near the phase boundary between the 
competing 6c-cycloidal and i?-type phases. 
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